Cold N+NH Collisions in a Magnetic Trap 
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We present an experimental and theoretical study of atom-molecule collisions in a mixture of cold, 
trapped atomic nitrogen and NH molecules at a temperature of ~ 600 mK. We measure a small 
N+NH trap loss rate coefficient of fc 1( ^~ N = 8(4) x 10~ 13 cm 3 s _1 . Accurate quantum scattering 
calculations based on ah initio interaction potentials are in agreement with experiment and indicate 
the magnetic dipole interaction to be the dominant loss mechanism. Our theory further indicates 
the ratio of N+NH elastic to inelastic collisions remains large (> 100) into the mK regime. 
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The study of low temperature molecular collisions and 
interactions is a rapidly expanding area of research at 
the interface of chemistry and physics [H, Experi- 
mental and theoretical techniques have been developed 
for studying a wide variety of phenomena, ranging from 
quantum threshold scattering 0, Q], inelastic atom - 
molecule collisions 0-0] , external field control of dipo- 
lar interactions in cold and ultracold molecules [1,13, 
to chemistry at cold 1(3- HJ and ultracold temperatures 
@j EH- The immense diversity of molecular structure 
and interactions is the cornerstone of many applications 
of cold molecules in quantum information science [lij ]. 
condensed-matter physics [Hj], precision measurement 



16l Il7| , cold controlled chemistry [18( , and astrophysics 
19j. Achievement of these applications relies on im- 
proved methods for cooling and further understanding 
of low temperature molecular collisions. 

Recently, several experimental techniques were devel- 
oped to study molecular collisions at cold and ultracold 
temperatures, beyond the reach of traditional molecular 
beam experiments [ioj . Ultracold ground state molecules 
with temperatures below 1 /xK can be created via co- 
herent state transfer from magneto-asssociated ultracold 
atomic gases, producing KRb [2l| and Cs2 22[. Al- 



ternatively direct cooling techniques such as buffer-gas 
cooling (23| and Stark deceleration Q provide more gen- 
eral methods for cooling a larger class of molecules, en- 
abling the production of cold molecules starting from 
room-temperature sources as an input, such as super- 
sonic beams Since the maximum molecular trap 
densities for these direct cooling techniques are, so far, 
10 8 — 10 9 enr 3 , the collision experiments with these 
molecular samples typically use a dense rare gas collision 
artner 0, Q or perform collisions at energies > 10 cm -1 
II > III- Further cooling and compression of molecular 
samples produced via direct cooling techniques is key to 
study cold molecule-molecule interactions and their ap- 



plications. One possible approach is to sympatheti cally 
cool a molecule via collisions with an atom 0, 124142 8j . 
Thus, studies of cold atom-molecule collisions not only 
uncover the underlying collision physics, but also may 
lead to important new methods for the production of a 
variety of ultracold molecules. 

A number of theoretical studies have focused on the 
optimal selection of atomic collision partners for sympa- 
thetic cooling of molecules [24 - 28 1 . Almost all of these 
studies focused on the alkali-metal atoms, which can be 
easily prepared at milli-Kelvin temperatures via laser 
cooling. However, accurate quantum scattering calcu- 
lations have shown that collision-induced inelastic relax- 
ation of molecules such as OH Q, ND 3 [25| and NH 
[2rj| with the alkali-metal atoms occurs rapidly, causing 
molecular trap loss and severely limiting the efficiency 
of sympathetic cooling. A recent theoretical study sug- 
gested that sympathetic cooling of NH molecules by col- 
lisions with laser-cooled alkaline earth atoms such as 
Mg might be possible at low magnetic field strengths 
271 ] . A related study indicated the possibility of sympa- 



thetic cooling of large molecules b y c ollisions with rare- 
gas atoms in an optical dipole trap [28[ . We have recently 
suggested that cold atomic nitrogen (N) can be used for 
sympathetic cooling of open-shell molecules, and demon- 
strated co-trapping of N atoms with NH molecules p3| . 
but the presence of 3 He buffer gas at 0.5 K precluded 
measurement of N+NH collisions in that experiment. 

In this Letter, we report the observation of cold (570 
mK) collisions between magnetically trapped, ground 
state open-shell atoms and ground-state polar molecules, 
N and NH, both species spin polarized. We measure 
a small N+NH trap loss rate coefficient of ^| ( ^^ NH ' ) = 
8(4) x 10~ 13 cm 3 s _1 . To interpret experimental obser- 
vations, we carry out accurate ab initio and quantum 
scattering calculations of N+NH collisions in a magnetic 
field. Our theoretical results agree with the observations 
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FIG. 1: NH trap decay taken at a trap depth of 3.9 T and 
cell temperature 570 mK. The solid line is a fit to the form 
n N H = aexp[-r N H^]. 



FIG. 2: NH loss rate vs cotrapped N density. The data 
are fitted to Eq. (1), yielding a value of fc 1( ^j~ NH ' = 8(4) x 
10" 13 cmW 



and predict high clastic-to-inelastic ratios (7) over a wide 
range of temperatures (from IK down to about 1 mK), 
suggesting the feasibility of sympathetic cooling of NH 
molecules by collisions with N atoms in a magnetic trap. 
Recent unpublished theory work by another group has 
made similar theoretical claims (29| . 

The apparatus used to co-trap atomic nitrogen and 
NH is similar to that described in our previous work 23, 
[30l | . A pair of super-conducting solenoids produce a 4 T 
deep spherical quadrupole magnetic trap. In the bore of 
the solenoids resides a cryogenic buffer gas cell held at 
500 mK. N and NH are produced in a molecular beam 
using a DC glow discharge and enter the trapping region 
through a 1 cm diameter aperture in the buffer gas cell. 
N and NH thermalize with 3 Hc buffer gas to 500 mK and 
fall into the trap. The initial buffer gas density at time 
of loading of 10 15 cnr 3 is produced using a cryogenic 
reservoir with fast actuating valve pSCjj to inject buffer 
gas through a 3.8 cm diameter aperture in the cell. After 
trap loading, the valve is closed, and the buffer gas exits 
primarily back out through the larger aperture, yielding 
background helium densities of 10 12 cm -3 , corresponding 
to NH trap lifetimes of several seconds. 

Detection of trapped N is performed using two-photon 
absorption laser induced fluorescence (TALIF) from the 
ground (2p 3 ) 4 S 3 / 2 state to the excited (3p) 4 S 3 / 2 state at 
96750 cm -1 . NH is detected via laser induced fluores- 
cence, as described in [23j]. We initially load NH with 
densities on the order of 10 8 cm~ 3 . Accurate knowl- 
edge of the absolute NH density is unnecessary as the 
NH density is much smaller than the N density in all our 
measurements. To determine the N density, the TALIF 
signal is calibrated using N+N collisional loss measure- 
ments [3l|. We load 5 x 10 11 N atoms into the trap with 
peak densities of 10 12 cm" 3 . 

Figure [T] shows a typical NH trapping decay. At time 
t = s, N atoms and NH molecules are co- loaded into the 
magnetic trap. For the first 2 s the NH trap loss is rapid 
due to the collisions with background helium gas during 
the pump out of buffer gas from the cell. By t = 5 s 



the NH trap loss reaches a steady rate, with typical 1/e 
lifetimes of about 3 s. Since the excitation laser for the 
atomic N TALIF detection causes additional loss of the 
trapped NH, we measure the trapped atomic N density 
at t = 15 s. 

To measure N+NH collisions, we observe the NH trap 
loss over a range of co-trapped N densities. The NH 
trap loss is fit between t = 5 and 10 s to an exponential 
decay to extract the total NH trap loss rate, Fnh- The 
cotrapped N density is varied by changing the ratio of 
molecular beam process gases N 2 and H 2 between (90%, 
10%) to (3%, 97%). Figure H shows the total NH loss 
rate, Fnh, versus cotrapped nitrogen density. The solid 
line in Fig. [5] is a fit to the equation 
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where tin is the peak nitrogen density, Tue is the NH 
loss rate attributable to collisions with background he- 



lium gas, and k 



.(N+NH) 



is the N+NH loss rate coefficient. 



The factor of 14 arises from averaging the N and NH 
densities over the volume of the magnetic trap. We find 
from this fit that k^ NH) = 8(4) x 10~ 13 cir^s" 1 . The 

uncertainty in fcf c ^ NH ' 1 is dominated by the uncertainty 
in our nitrogen density calibrations [3l| . The N+NH loss 
rate coefficient has contributions from both elastic (evap- 
orative) and inelastic N+NH collisions [3l|, but their in- 
dividual contributions cannot be determined by a single 
measurement. Therefore, our measurement of fc^^ NH '' 
is an upper limit on the inelastic N+NH rate coefficient 
and could be lower by up to a factor of 50%. 

To interpret the experimental observations and explore 
the possibility of sympathetic cooling of NH molecules 
by collisions with co-trapped N atoms, we performed rig- 
orous quantum scattering calculations of inelastic relax- 
ation in N+NH collisions in an external magnetic field. 
The Hamiltonian of the atom-molecule collision complex 
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FIG. 3: Contour plot of the ab initio interaction PES for the 
fully spin-polarized S = 5/2 state of N+NH calculated in this 
work. The NH bond distance is fixed at its equilibrium value 
r = 1.958 do. Energies are in units of cm - . 



may be written (h = 1) 
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where fi is the reduced mass of the complex, R = \R\ is 
the N+NH separation vector, r is the internuclear dis- 
tance in NH, 6 is the angle between the vectors R and 
r, Vs(R,r,8) is the interaction potential of the N+NH 
collision complex, S is the total spin of the complex, and 
Ms is its projection on the magnetic field axis. The last 
two terms in Eq. ([2j describe non-interacting collision 
partners in the presence of an external magnetic field of 



strength B |3j 



The term proportional to R 



the magnetic dipole interaction [31 



represents 



The interaction of NH( 3 £) molecules with N( 4 S 3 / 2 ) 
atoms gives rise to three adiabatic potential energy sur- 
faces (PESs) with S = 1/2, 3/2, and 5/2. Since in 
our experiments both atoms and molecules are confined 
in a permanent magnetic trap, the incident collision 
channel is the maximally spin-stretched Zeeman state 
|M Snh = 1) ® \M Sn = 3/2) with S = 5/2. We neglect 
the weak couplings between the electronic states of dif- 
ferent S arising from the fine-structure terms in i?NH 
and off-dia gonal m atrix elements of the magnetic dipole 
interaction 
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To evaluate the PES for the S = 5/2 electronic state 
of N-NH, we use the state-of-the-art partially spin re- 
stricted coupled cluster method with single, double, and 
non-iterative triple excitations [RCCSD(T)] [33| using 



quadru ple- zeta basis set (aug-cc-pvqz) of Dunning et 
al. [34 - 36 1 augmented with 3s3p2e?2/lg bond functions 
placed in the middle of the intermolecular distance. A 
contour plot of the calculated PES in shown in Fig. 3. 
The PES has a global minimum 87.83 cm -1 deep in the 
linear N-HN configuration (R = 7.02ao, = 0) and a 
secondary minimum of 77.52 cm -1 in the HN-N configu- 
ration (R = 6.61ao, = 180°) separated by a barrier 39.6 
enr 1 high located at 9 = 92°. Compared to the PES for 
the He-NH interaction [13|, our calculated N+NH PES 
is both deeper and more anisotropic. 

To solve the scattering problem, we expand the wave 
function of the collision complex in a direct-product ba- 
sis set \NM N }\S N nMsn H )\S N Ms N )\lm e }, where N is the 
rotational angular momentum of NH, and M/v, Ms Na , 
Ms N , and mi are the projections of N, SnHj Sh, and £ 
on the magnetic field axis. This expansion results in a 
system of close-coupled (CC) equations parametrized by 
the matrix elements of the Hamiltonian @, which can 
be readily evaluated analytically given the direct-product 
structure of the basis [32|. We solve the CC equations 
numerically for each value of the total angular momen- 
tum projection M = AIn + A/s NH + Ms N + m e f° r colli- 
sion energies between 10 and 1 cm -1 and extract the 
S'-matrix elements and scattering cross sections. Large 
basis sets with N = — 5 and £ = — 8 are used to ensure 
that the results arc converged to < 5%, resulting in 2906 
CC equations for M = 0. 

Figure 4(a) shows the cross sections for elastic energy 
transfer and inelastic relaxation in N+NH collisions. The 
cross sections increase with decreasing collision energy 
before reaching a maximum at Eq ~ 2 mK, which we 
identify as a shape resonance supported by the centrifu- 
gal barrier with £ = 1 in the incident collision channel. 
At Eq < 0.1 mK the inelastic cross sections assume the 
characteristic l/Ec dependence on collision energy, and 
the elastic cross sections become constant, according to 
the Wigncr threshold law for s-wave scattering. By ther- 
mally averaging the cross sections in Fig. 4(a), we obtain 
elastic and inelastic rate constants at the experimental 
temperature of 570 mK of fc cl = 2.2 x 10~ 10 cm 3 s _1 , 
k in = 4.1 x 10~ 13 cm 3 s _1 . The fraction of elastic N-NH 
collisions that lead to trap loss is on the order of 3 x 10~ 3 , 
leading to a total trap loss rate coefficient fe™j^!S 



loss — theory 

enr's^ 1 , in good agreement with the experi 
mental value of fcf* +N *2 - '« + ^ - 1 n " 13 — 3 °" 1 



11 x 10" 13 — 3 = 



loss-experiment = (8±4)x 10 13 CI 

The ratio 7 = k e i/fa n is shown in Fig. 4(b) as a function 
of temperature. The ratio remains large (7 > 100) over 
the temperature range ~10 mK - IK, which indicates 
that NH molecules can be sympathetically cooled by elas- 
tic collisions with spin-polarized N atoms down to the 
milli-Kelvin regime. As shown in Fig. 4(b), an increasing 
magnetic field suppresses inelastic relaxation, so that 7 
remains high (~50) even at T = 1 mK. This observation 
indicates that applying a strong uniform magnetic field of 
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FIG. 4: (a) Cross sections for elastic scattering (circles) and 
inelastic relaxation (triangles) in NH + N collisions calculated 
as functions of collision energy at a magnetic field of 0.1 T. 
Also shown is the cross section calculated with the magnetic 
dipole interaction omitted (dashed line), (b) Thermally av- 
eraged ratios of the rate constants for elastic scattering and 
inelastic relaxation as functions of temperature for B = 0.1 
T (triangles) and B = 1 T (diamonds). 



order 1 T may be used to stabilize spin-polarized atom- 
molecule mixtures against collisional losses, thereby en- 
hancing the efficiency of sympathetic cooling. 

Inelastic collisions of 3 S molecules with open-shell 
atoms like N can occur due to (i) indirect couplings in- 
duced by the anisotropy of the atom-molecule interac- 
tion potential and the spin-spin interaction between the 
N = and N = 2 rotational states [H, [32j , and (ii) direct 
couplings between atomic and molecular Zeeman levels 
induced by the long-range magnetic dipole interaction^ • 
The magnetic dipole interaction couples the Zeeman lev- 
els directly 31 1 , and can thus be more efficient in induc- 
ing inelastic relaxation than the indirect couplings. Fig- 
ure 4(a) shows that omitting the magnetic dipole inter- 
action from scattering calculations reduces the inelastic 
cross sections by a factor of ^10, confirming that Zeeman 
transitions in N+NH collisions are indeed driven by the 
magnetic dipole interaction. In our experiments, both N 
and NH collision partners are fully spin-polarized, so the 
chemical reaction N + NH — >• N2 + H is spin-forbidden 
and can only proceed via non-adiabatic transitions be- 
tween different electronic states of the N-NH complex 
in the entrance reaction channel mediated by the fine- 
structure and magnetic dipole couplings (sec Eq. 2). Our 



observed value for k 



(N+NH) 



is much smaller than the 



measurements are consistent with theoretical predictions 
which do not account for the reaction channel, we con- 
clude that chemical exchange processes in spin-polarized 
N+NH mixtures occur at a slow rate, and do not con- 
tribute to the observed trap loss dynamics. This im- 
portant finding shows that inelastic relaxation in N+NH 
collisions occurs via the same mechanism as dipolar re- 
laxation in spin-polarized atomic gases (3l| . 

In conclusion, we have measured a small N+NH trap 
loss rate coefficient of fc|^ NH ^ = 8(4) x 10~ 13 cm 3 s _1 at 
a temperature of ~ 600 mK. To interpret experimen- 
tal observations, we have carried out accurate ab initio 
quantum scattering calculations of Zeeman relaxation in 
N+NH collisions in a magnetic field and find theory and 
experiment to agree. Our calculations show that the ra- 
tio of N+NH elastic to inelastic collisions remains large 
(> 100) over the temperature range ~10 mK - IK, which 
indicates that it may be possible to sympathetically cool 
NH molecules down to the milli-Kelvin regime via elas- 
tic collisions with spin-polarized N atoms. It remains to 
be seen whether this conclusion holds for other paramag- 
netic molecules of interest such as the highly polar CaH 
or SrF in their electronic ground states of 2 E symmetry. 
If it does, it may be possible to create large samples of 
these molecules via collisional cooling with N atoms in a 
magnetic trap. 
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calculated reaction rate for spin-unpolarized reactants 
(Reaction ~3x 10" 11 cm 3 /s at T = 1 K) [38|. Since our 
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